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. . ABSTRACT 
>>. 

■ The Kerr-lens mode locking ability and the ultrashort pulse characteristics are analyzed numerically for the Cr-doped 
ZnTe, ZnSe, ZnS active media. The advantages of these materials for the femtosecond lasing within 2-3 fim spectral 
range are demonstrated. 
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Compact diode-pumped sources of the femtosecond pulses tunable within the wavelength range between 2 and 3 
/im are of interest for various applications, such as laser surgery, remote sensing and monitoring, spectroscopy of 
^ semiconductors etc. To date only cryogenically operated Pb-salt diode lasers, optical parametrical oscillators and 
^ difference-frequency convertors were available for the operation in this spectral range. Th erefo re the possibility of 
' ^ ' the direct mid-IR lasing from the new class of the transition- metal doped chalcogenides has attracted much 

J>^, attention. The impressive advantages of these media are room-temperature operation between 2 and 3 /j,m, possibility 
'■ of direct diode pumping, high emission and absorption cross-sections, negligibly low excited-state absorption and, as 
consequence, low thermal load (the basic laser material characteristics will be described in the next section). The 
most remarkable examples of such lasers are Cr^+-doped ZnSe, ZnS and ZnTe. To date the following achievements 
Pn| ' for these media have been demonstrated: 1) for Cr:ZnSe CW operation with over 1.7 W power P over 1100 nm 
'i^ tunability P over 350 nm tunable diode-pumped CW operation ^ a ctive mode locking ^ and active modulator 
] assisted passive mode lockingl^were achieved; 2) for Cr:ZnS pulsed I^ESl and tunable CW operation El were obtained. 
^"^^ Cr:ZnTe, which is a member of the considered media class, remains unexplored. 



(N 
O 



In spite of the numerous advantages, there exist some obstacles for femtosecond pulse generation from these lasers. 
As they are the semiconductors, i.e. possess a comparatively narrow band-gap, the nonlinear refraction in the active 
crystal is large (see below). Hence the self- focusing has a comparatively low threshold which in the combination with 
the self -pha se modulation produce a tendency to the multiple pulse operation in the Kerr-lens mode locking (KLM) 
Q ' regime Thus, there is a need in the study of the KLM stability limits and methods of the stability enhancement 
' 55 ' considered lasers. Moreover, the large absorption and emission cross sections tend to the stabilization of the 

J>-^' CW operation, which prevents the KLM self-start. 



In this paper we present an analysis of the nonlinear refraction in the Zinc-chalcogenides (by the example of 
ZnSe). In the combination with the lasing properties this defines the main requirements to the KLM optimization. 
Then, the results of the numerical optimization of the KLM aimed at the multipulsing suppression and taking into 
^> ' account a strong saturation of the Kerr-lens induced fast absorber are presented. We demonstrate the possibility of 
^ the few optical cycle pulse generation from the Cr-doped Zinc-chalcogenodes. The problem of the self-starting KLM 
5^ is analyzed on the basis of the generalized momentum method taking into account real-world laser configurations, 
gain guiding, soft-aperture, thermo-lensing and other properties of the considered laser materials. The presented 
models are quite general and can be applied to the overall optimization of the different KLM lasers. 
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2. DISTRIBUTED MODEL AND BASIC PARAMETERS 



Simulation of the KLM can be based on the two quite different approaches. First one supposes the full-dimensional 
modelling taking into account the details of the field propagation in the laser cavity The minimal dimension 
of such models is 2+1 and they allow the description of the spatio-temporal dynamics of the ultrashort pulses and 
their mode pattern. Its main disadvantages are a large number of the parameters resulting in ambiguity of the 
optimization procedure and complexity of the interpretation of the obtained results. Second approach is based on 
1+1 dimensional model in the framework of the so-called nonlinear Ginzburg-Landau equation which describes 
the KLM as an action of the fast saturable absorber governed by the few physically meaningful parameters, viz., 
its modulation depth 7 and the inverse saturation intensity a. This method allows the analytical realization in the 
week-nonlinear limit however in the general case the numerical simulations are necessary. We shall base our 
analysis on the both approaches. 

In the beginning let us consider the master equation describing the ultrashort pulse generation in the KLM 
solid-state laser: 



da{z, t) 
dz 




i5[\a{z,t)\^ - —^^\a{z,t)\^ 



a\a{z,t)\^ V ' ^odt 



a{z,t), 
(1) 

where a{z,t) is the field amplitude (so that |ap has a dimension of the intensity), z is the longitudinal coordinate 
normalized to the cavity length (thus, as a matter of fact, this is the cavity round-trip number), t is the local time, a is 
the saturated gain coefhcient, p is the linear net-loss coefficient taking into account the intracavity and output losses, 
tf is the group delay caused by the spectral filtering within the cavity, are the m-order group-delay dispersion 
(GDD) coefficients, S = lgn2i^o/c = 2TTn2l g / {Xqu) is the self-phase modulation (SPM) coefficient, ujq and Aq are the 
frequency and wavelength corresponding to the gain band maximum, n and n2 are the linear and nonlinear refraction 
indexes, respectively, Ig is the double length of the gain medium (we suppose that the gain medium gives a main 
contribution to the SPM). The last term in Eq. describes the self-steepening effect and for the simplification 
will be not taken into account in the simulations. As an additional simplification we neglect the stimulated Raman 
scattering in the active medium ^ These two factors will be considered hereafter. 

The gain coefficient obeys the following equation: 

daiz,t) Ip \a{z,t)f a{z,t) 

= a. («^.. - a{z, t)) — - agaiz, i)^- —. (2) 

Here ctq and ag are the absorption and emission cross-sections of the active medium, respectively, Tr is the gain 
relaxation time, Ip is the absorbed pump intensity, tOp is the pump frequency, a^ax = o'gNglg is the maximum gain 
coefficient, Ng is the concentration of the active centers. The assumption Tp <^ Tcav {Tp is the pulse duration, Tcav 
is the cavity period) allows the integration of Eq. (0) . Then for the steady-state gain coefficient we have: 



where Eg = hujp/crg is the gain saturation energy fiux, E = \o,{t)\^dt is the pulse energy. 

For the numerical simulations in the framework of the distributed model it is convenient to normalize the time 
and the intensity to tf = Ag/ (AAc) and 1/S, respectively (AA is the gain bandwidth). The simulation were performed 
on the 2^^ x 10^ mesh. Only steady-state pulses were considered. As the criterion of the steady-state operation we 
chose the peak intensity change less than 1% over last 1000 cavity transits. 

The KLM in the considered model 1^131 is governed by the only four basic parameters: a — p, f32, J, and a. This 
allows unambiguous multiparametric optimization. In the presence of the higher-order dispersions, the additional 
(3rn parameters appear. This complicates the optimization procedure, but keeps its physical clarity. 

Now let us consider the basic active media parameters. 



Table 1. Material parameters of the Cr-doped Zinc-chalcogenides. 



Medium 


Ao, A*ni 


AA, nm 


Aq, ^m 


(Ta, 10 i''^ cm-^ 


(7g, 10 cm'' 


n 


712, 10 esu 


Tr, ^s 


Cr:ZnSe 


2.5 


880 


1.61 


8.7 


9 


2.44 


29 - 92 


6-8 


Cr:ZnS 


2.35 


800 


1.61 


5.2 


7.5 


2.3 


16 - 30 


4-11 


Cr:ZnTc 


2.6 


800 


1.61 


12 


20 


2.71 


80 - 150 


3 



3 




Figure 1. Calculated wavelength dependence of n2 for ZnSe. Solid and dashed curves present the maximum and 
minimum estimations based on Eq. Points give the experimental data for polycrystalline samples, squares - for 
monocrystals with different orientations, crosses - unknown crystalline structures. 

It should be noted that the n2 values were experimentally measured at the wavelengthes which are smaller than 
Ao (see Fig. Therefore we have to estimate theirs at the generation wavelength. Such estimation can be obtained 
from the two ™ or four bands models of semiconductor. The most simple estimation results from the former 
model, which gives the following formula: 

"2 (esu) = K — , (4) 

where K = (0.5 - 1.5)xl0~® and Ep = 21 eV are the material independent constants. Eg is the band-gap width in 
eV, G is the form-factor. Using for K the values 0.86-1.5 xlO~* and for Ep the values 21-24 eV, we obtained the 
estimations presented in Tabled The calculated n2 dispersion for ZnSe is shown in Fig. ^P^l 

We note that the semiconductor nature of these active media results in the extremely high nonlinear refraction 
coefficients in the comparison with Ti:sapphire (1.2x10^^^ esu). Also, we have to note that there exists the orien- 
tational dependence of 712, which can be especially large for the birefringent ZnS. But even for the cubic ZnSe, the 
orientational variation of rt2 can reach 70% (see Fig. 

In addition to the large n2 the considered media possess the strong quadratic-nonlinearity due to absence of the 
center of inversion (for ZnSe and ZnTe). For example, dag = 33 pm/V for ZnSe. However, as a result of the cubic 
structure, there exists the strong orientational dependence of the effective second-order nonlinearity (see Fig. [JJ, 
so that the maximum effective nonlinear coefficient is only 0.17c?36. The coherent length for this media is only 23 
/im that strongly reduces an efficiency of the second-harmonic generation. Nevertheless, the cascaded second-order 




Figure 2. (a): Calculated oricntational dependence of the third-order rt2 for ZnSe 1^ scaled to rij*"^^' and (&): 
orientational dependence of the second-order nonlinear coefficient scaled to its maximum value. 
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Figure 3. Calculated wavelength dependence of the cascaded n2 for ZnSe. 

process can contribute to the effective 712 .I2E11 Our calculation demonstrate (see Fig. that this contribution 

is small in the comparison with the third-order 712 . 

As a result, the simulation parameters corresponding to the above introduced normalizations are summarized 
in Table 12 P — aaTcavIp/ {fi^a), e = ^f^^^/^'j Pcr—'^nQ / (n2k'^) and Isat—hv^liugTr) are the critical power of 
the self-focusing (beyond the parabolical approximation for the Gaussian beam) and the gain saturation intensity, 



respectively. 



Table 2. Simulation parameters. Ig— 2x0.3 cm, — ns, 2 W pump power, 100x100 fim^ pump mode. We used 
the average value of n2 and the minimum of Tj.. 



Medium 


^max 


712, 10-^" cmVW 


Per, MW 


S, IQ-''^ cmVW 


e, 10-* 


Isat, kW/cm^ 


t/, fs 


P, 10"^ 


Cr:ZnSe 


5 


105 


0.75 


87 


4.9 


15 


3.8 


1.1 


Cr:ZnS 


5 


42 


1.5 


32 


10 


28 


3.7 


0.62 


Cr:ZnTe 


5 


178 


0.52 


314 


3.8 


13 


4.5 


1.9 



3. DISTRIBUTED MODEL: RESULTS AND DISCUSSION 

First of all we have to consider the meaning of the optimization procedure. The numerical simulations demonstrate 
that there exist some values of the saturation parameter tr, which provide the near chirp-free pulse generation. These 
values of a can be considered as "optimal". However, there is an additional factor, which has to be taken into 
account: the pulse shortening is possible by the a growth and P2 ^0. We consider the values of (32 and a as optimal 
if they correspond to the generation of the shortest pulse. The main obstacle on the way of the pulse shortening 
is the multiple pulse generation.^ The strong tendency to the multipulsing results from 1) decrease, 2) P 
increase, 3) fast absorber saturation favored by the large a. The relatively large values oiag and Aq for the considered 
media result in the increase of the first parameter. However, the large absorption cross sections and SPM coefficients 
increase two later parameters. If the growth of the first parameter corresponds to the gain saturation and thereby to 
the stabilization of the pulse against the laser continuum growth, the larger P and a initiate the rise of the continuum 
and the excitation of the perturbations inside the ultrashort pulseP^l Hence, the optimization means the multiple 
pulses suppression allowing the pulse shortening. 

Additional limitations on the pulse shortening result from the achievable values of the saturation parameter a. 
In the KLM lasers this parameters is governed by the cavity alignment: the shift towards the cavity stability zone 
increases a (see next section). So, the highest values of a are reached in the immediate vicinity of the cavity stability 
boundary. This requires a too thorough cavity optimization and can not be considered as operational. Hence, the 
optimization aimed at the pulse shortening is based on the variation of all four parameters of master Eq. JQ) and is 
constrained by the above described reasons. 

Fig-Elshows the parameters of Eq. Q|, which for the fixed a — p, 7 and a correspond to the minimum achievable 
/32, where the pulse width has a minimum. The further P2 decrease results in the multipulse operation. Thus, the 
points in Fig. 01 lie on the boundary of the stable single pulse operation. There is a set of the general tendencies 
characterizing this boundary. 

There exists a limited on a ~ p range of the stable single pulse operation, which expands as a result of the a and 
7 growth. The increase of a shifts this region towards the smaller a — p. However, when a >100 the transition to 
multipulsing is possible. As every point in Fig. 0] corresponds to the fixed dimensionless E (see Fig. EJ, the choice 
of the appropriate level of a — p for the fixed laser configuration (i.e. fixed 7, a and (32) is realized by the variation 
of P (pump, see Eq. ^) as well as p (output loss). It is possible also to change a by the change of amax due to 
variation of the crystal length or the active ions concentration or by the change of due to variation of the active 
ions concentration. Note also, that the Ig decrease increases e (due to the 6 decrease), which describes the "strength" 
of the gain saturation relatively the SPM. Hence, e increase expands the stability region towards the higher pump 
and allows the higher pulse energies (because they oc 1/S). 

The E decrease for the fixed a — p, which takes a place for (32 — >0 (Fig.^l, has to be accompanied by P, amax, 
Tr decrease or by e, p increase in order to prevent from the continuum amplification. The last is the main source 
of the pulse destabilization and suppresses the single pulse operation in the vicinity of zero GDD. Hence, the pulse 
generation for a — p — 7 >0 is not possible, a — p — 7 =0 corresponds to the specific hybrid regime with the coexistent 
pulse and continuum 

Since the approach of the GDD to zero has to be accompanied by the pump decrease, this results in the growth 
of the a required for the pulse stabilization (Fig. This can demand too thorough cavity alignment. Moreover, 




Figure 6. Regions of the a and /?2 parameters providing tlie shortest pulses. 

for the large a we need the larger minimum \[32\ providing the single pulse operation so that the dependence of the 
minimum |/32| on the a for the fixed a — p has a parabolic- like form P^l 

The most interesting features are the shift of the stability region towards the smaller cr (Fig. 0)) and the pulse 
shortening as a result of the 7 increase. For example, the minimum pulse duration Tp for 7 =0.03 is 7tf whereas for 
7 =0.05 it is 5tf (this is 19 fs for Cr:ZnSe and Cr:ZnS and 23 fs for Cr:ZnTe). The bad news here is the need for the 
hard-aperture KLM to provide the larger modulation depth. This reduces the KLM self-starting ability. 

The regions of the parameters allowing the shortest pulses are shown in Fig. El The a — p increase reduces the 
minimum a parameter producing the shortest pulses. However, the region of their existence shortens on cr. 




Figure 7. Pump threshold of KLM related to that for CW. Pulse duration is Tp =7 t/. a =0.071, p =0.05, 7 =0.03, 
ctrnax =5, Other parameters correspond to Table[21 



Let us consider the concrete parameters of Table |21 The pump thresholds allowing 7tf pulse durations are shown 
in Fig. 13 The threshold decreases from CriZnS through Cr:ZnSe to CriZnTe that results from the e decrease. This 



means that the SPM becomes stronger relatively the gain saturation under this transition. As a result, the tendency 
to the pulse destabilization intensifies and this demands to reduce the intracavity power by means of the pump 
decrease. 




Figure 8. Intracavity pulse energy fluxes. The parameters correspond to Fig.[7| 

Thus, the e decrease reducing KLM threshold turns in the intracavity pulse energy decrease (Fig. (SJ. The 
highest value of e for Cr:ZnS resulting from the large dg in the combination with comparatively small S produces the 
stabilization of the shortest pulses with the highest energies. 

Note, that the larger absorption cross-section for Cr:ZnTe results in the highest absorbed pump energy for the 
fixed pump intensity and mode area. This is a positive factor for the KLM threshold lowering. However, this can 
be a negative factor, when the SPM is the source of the pulse destabilization because the additional efforts for the 
intracavity power decrease are necessary (see also next section). 

At last, we consider the contribution of the third-order dispersion, which can be large for the lasers under 
consideration. There are the technological troubles in the use of the chirped-mirror technique for the dispersion 
compensation in the mid-IR due to high value of Aq. Therefore the usual technique utilizing the prisms for the 
dispersion control can be useful in the considered situation. As a result, the third-order net-dispersion coefficient 
1/33 1 increases. 

For the simulation we choose =-5900 fs^. As it can be seen from Fig. ^ the shape of the stability boundary 
does not change in the comparison to f]^ =0 .^^ However, the minimum pulse duration increases from 5tf for —0 
to 9tf (34 fs for Cr:ZnS, Cr:ZnSe and 40 fs for Cr:ZnTe). The additional effect is the pronounced (up to 140 nm) 
Stokes shift of the peak wavelength (Fig. |51shows this shift on the stabilit y boun dary for Cr:ZnSe and Cr:ZnS). This 
shift is typical also for such IR lasers as Cr:LiSGaF, CriLiSAF, Cr'^+ : YAG l ^^ l ^^ l and can reduce the pulse energy due 
worse overlap between gain band and pulse spectrum. However, for the media under consideration the wavelength 
shift is small in comparison with the full gain band width and the energy decrease is not critical. 

4. TIME-SPATIAL MODEL: SELF-START ABILITY OF THE KERR-LENS MODE 

LOCKING 

In the previous section we considered the ultrashort pulse stability on the basis of the distributed 1-fl-dimensional 
model. However, this model can not provide answers to the basic questions: what is the self-start ability of the 
Kerr-lens mode locking and what has the real-world laser configuration to be? 

These questions require an analysis on the basis of the time-spatial model. The simplest way is the assumption 
of the Gaussian spatial distribution for the laser beam reducing problem to the l-|-2-dimensions. The free-space 
propagation of the Gaussian beam can be considered on the basis of the usual ABCD-matrix formalism, while the 
propagation inside the nonlinear active medium is described by the following equation: 



y=0.05 




Figure 9. Stability boundary for the single pulse operation in the presence of the third-order dispersion and the 
Stokes spectral shift of the pulse spectrum on this boundary. 
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Here P2 and t'j, are the group-velocity dispersion and the inverse group-velocity delay coefficients (for ZnSe we used 
/3^2=2054 fsVcm and t'^^lS fs/cm). 

The left-hand side of Eq. (0 describes the non-dissipative factors: thermo-lensin^^H (^=fc^a.(^p^ exp (— (^2;)/ 
{4:TTnoKth), k is the wave number, ^S- is the coefhcient of the refractive index thermo-dependence (5.35x 10~^ii'~"'^ 
for ZnSe), ( is the loss coefficient at the pump wavelength, Pa is the pump power, Kth is the thermo-conductivity 
coefficient (0.172 W K~^cm~^ for ZnSe)); diffraction (in the cylindrically symmetrical case); group- velocity dispersion 
and self-phase modulation (providing self-focusing for radially varying beam, x=i^2k/no). The right-hand side of Eq. 

describes the dissipative factors inside the gain medium: radially varying gain (providing gain guiding and soft 
aperture action, a and Wp are the saturated gain coefhcient and the pump beam size, respectively); spectral filtering 
caused by the gain band profile. 

It is convenient to rewrite Eq. ^ in the following way: 



P T 

(6) 



where v=Et:w'^ / {2PgTp) {Pg is the generation power), w is the generation mode beam size, v—^'^j^ for the pulse 
with the Gaussian time-profile, 2 for the sec/i-shaped pulse and 1 for the CW (in the latter case Tp=Tcav)- 

The crucial simplification in the analysis of Eq. ij^l is based on the so-called aberrationless approximation: the 
propagating field has the invariable spatial-time profile, which is described by the set of the z-dependent parameters. 
In the non-dissipative case this approximation allows the variation al appro ach providing rigorous description of the 
Gaussian beam propagation outside the parabolical approximation P'^ l ^'^ l ^^ l 

In the dissipative case we have to use the momentum method ."^ However, in contrast to Refs. EmU we shall 
consider the momentums resulting from the symmetries of Eq. (|5|). The a — > aexp(i(/)) invariance, the transverse 
and time translating invariance suggest the following momentums 1221 




oo 



Like the variational approach we can substitute to Eqs. I0J the trial solution describing the ultrashort pulse. If 

we take the Gaussian time-spatial profile a (z, r, t) = W{r) exp (^G{r) — ^vJlr)'^ + ib(r)r^ — -^^pp + i'4}{r)t^^ {W{r) is 

the complex amplitude, 2w''^=w^ , G{r) is the pulse amplification parameter excepting the geometrical magnification 
for the Gaussian beam), the equations describing the evolution of the pulse and beam parameters arj^ 
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2t) 

7{zf 



dh{z) _ 2^ ^ 2h {zf ^ y/2PQe^^^''^ 



dip (z) 
dz 



dz 



I 



r{z) 



k 

4 



nPcrkw (z)'^ 2fcw' (z)^ 



2\ 



(8) 



r(z)^ -K PcrW {z)'"^ T {z)^ 
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where Pq and w'^ are the power and the beam size before the active medium, respectively. 

Eq. © in the combination with the ABCD-formalism allows defining the stability regions for CW, single- and 
multiple pulsing (the latter requires the trivial generalization of the trial function). Such regions can predict the 
real-world laser configurations providing the self-starting Kcrr-lcns mode locking. 

As it is seen from Table 2 the main features of the considered active media are the comparatively low Per (e.g. 
Pcr=2.6 MW for Tiisapphire) and the extremely low Isat (300 kW for Ti:sapphire). This causes the stabilization 
of the CW against the mode locking due to the strong gain saturation. However, as it is seen from Eqs. (jHl El 
(note that the gain saturation is intensity-dependent, but the self-focusing is power-dependent), we can enhance the 
tendency to the mode locking due to the suppression of the gain saturation. This can be achieved by the growth of 
the generation mode within the active medium owing to, for example, the cavity shortening or the decrease of the 
folding mirrors curvature (Fig. I10|l . At the same time we have to keep the moderate level of the absorbed pump 
power and to use the asymmetrical cavity in order to separate the stability zones. The preferable configurations 
correspond to the pulse operation without the GW. This occurs between the GW stability zones (see Fig. I1(J|I . 



5. CONCLUSION 

In conclusion, we presented the models, which can be usable for the numerical optimization of the KLM lasers. 
On the basis of these models, the KLM abilities of the Gr-doped Zinc-chalcogenides were estimated. It was shown, 
that the strong SPM inherent to these media and destabilizing the single pulse operation can be overcomed by the 
choice of the appropriate GDD, pump, modulation depth and saturation parameter of the Kerr-lensing induced fast 



saturable absorb(^r. As a result, the Cr:ZnTc possesses the lowest KLM threshold, however strong SPM eonstrains the 
achievable pulse power for this laser. The best stability for the highest energy flux and the shortest pulse duration (19 
fs) are achievable in Cr:ZnS. Cr:ZnSe. The presence of the third-order dispersion increases the minimum achievable 
pulse durations up to 34 - 40 fs and causes the strong (up to 140 nm) Stokes shift of the generation wavelength. 
However, the latter effect does not reduce noticeably the pulse energy. The Kerr-lens mode locking self-start ability of 
the considered active media is reduced by the strong gain saturation so that the cavity tuning providing the ultrashort 
pulsing differs essentially from that for the Ti:sapphire. We have to avoid over-pumping and over-shortening of the 
generation mode. As a result, the short and asymmetrical cavities with the comparatively large curvature of the 
folding mirrors are preferable. On the whole, the Cr-doped Zinc-chalcogenides have the prospects for the sub-50 fs 
generation that amounts to only few optical cycles around 2.5 /um. 
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Figure 10. Stability ranges for the CW (black filled circles) and mode locking (open green circles), a and c are the 
cavity arms, h is the folding distance, h\ is the distance of the active medium facet from the curved mirror forming 
the shorter cavity arm, / is the focal length of the curved mirrors. p=0.05. 



